Phenomenology of lepton flavor violation in 
2HDM(III) from (g — 2) and leptonic decays. 

Rodolfo A. Diaz, R. Martinez, and J. -Alexis Rodriguez 
Universidad Nacional de Colombia. 



O 

o 

(N 
o . 

<D . Departamento de Fisica. Bogota, Colombia 



February 1, 2008 



' Abstract 

> ' 

' We constrain some lepton flavor violating vertices in the context of 

the two Higgs doublet model type III, assuming that the lightest scalar 
Higgs mass m h o is about 115 GeV. Specifically, based on the g — 2 muon 
QO ' factor and the decay width of /i — » cy, the following quite general bounds 

are obtained: 7.62 x 10" 4 < < 4.44 x 10" 2 , £ 2 T < 2.77 x 1(T 14 . 
CNl ■ Additionally, based on the processes r — * fiy, and r — > nnn, bounds on 

' £ TT and £ MM are also gotten, such constraints on these parameters still 

pH , give enough room for either a strong suppression or strong enhancement 

Q 1 on the coupling of any Higgs boson to a pair of tau leptons or a pair 

I ' of muons. Furthermore, upper limits on the decay widths of the leptonic 

decays r — > e-y, and r — > eee are calculated, finding them to stay far from 
the reach of near future experiments. 

> : 

■ 1 Introduction 

It is well known that processes involving Flavor Changing Neutral Currents 
(FCNC), are severely suppressed by experimental data, despite they seem not 
to violate any fundamental law of nature. On the other hand, Standard Model 
(SM) issues are compatible with experimental constraints on FCNC so far, with 
the remarkable exception of neutrino oscillations. Preliminary measurements of 
the flux of solar neutrinos Q yielded a result considerably smaller than the one 
expected from the standard solar model, from which the possibility of oscillation 
of neutrinos arose as a possible solution to the neutrino deficit observed Q. 
Later on, neutrino oscillations in matter were proposed to explain the deficit 
confirmed by SuperKamiokande ||]. Additional experiments with solar Q and 
atmospheric neutrinos || have supported the idea of the oscillations. 

On the other hand, this phenomenon implies the existence of mass terms for 
neutrinos, and the existence of family lepton flavor violation. Both implications 
lead us in turn to consider the existence of physics beyond the SM, since in the 
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SM neutrinos are massless and lepton flavor violating mechanisms are absent. 
Therefore, the increasing evidence on neutrino oscillations motivates the study 
of models with lepton flavor violation. There are some interesting scenarios in 
which FCNC and therefore LFV appears naturally. For example, the introduc- 
tion of new representations of fermions different from doublets produce FCNC 
by means of the Z-coupling S . On the other hand, they could arise at the tree 
level in the Yukawa sector if a second doublet is added to the SM II . Another 
interesting scenario appears in SUSY theories with R-parity broken jn^, be- 
cause FCNC coming from R-parity violation generates massive neutrinos and 
neutrino oscillations Jll| . FCNC in SUSY theories in the charged lepton sector 
have been studied as well. For instance, the decays /i — > ej and \i — ► 3e with 
polarized muons have been examined in the context of supersymmetric grand 
unified theories Q . 

Moreover, neutrino oscillations have inspired the study of Lepton Flavor 
Violating (LFV) processes in the neutral leptonic sector. For example, LFV 
processes in SU(5) SUSY models with right handed neutrinos have been exam- 
ined based on recent results of neutrino oscillation experiments On the 
other hand, ref Jl6| explores the generation of LFV by using a multi-Higgs dou- 
blet model with additional right handed neutrinos for each lepton generation, 
finding a non-decoupling behavior of some LFV amplitudes respect to the Right 
handed neutrino masses. 

Additionally, experimental search has imposed some upper limits for the 
branching of many processes involving LFV as in the case of — > /i + e~ fll7|] , 
K° L -> 7rV + e _ @, K+ -> Tr+^+e" @, n+ -> e+ 7 |§, [i+ -> e+e+er fl| 
and fi~N — ► e~N [^2|. As for experimental prospects, muon colliders offer a 
very appealing alternative because of its potential to test the mixing between 
the second and third (first) generations by processes like /j,fi — > jUr(er) mediated 
by Higgs exchange §, |§,|§. 

Perhaps the simplest framework to look for these rare processes is the Two 
Higgs Doublet Model (2HDM), which consists of considering the possibility of 
duplication of the SM Higgs doublet. Owing to the addition of the second 
doublet, the Yukawa interactions lead naturally to tree level FCNC unless we 
make additional assumptions. The most common one is the imposition of a 
discrete symmetry |l3| which in turn generates two independent sets of Yukawa 
interactions known as models type I and II. Notwithstanding, we can still have 
the possibility of preserving all terms in the Yukawa sector (known as model 
Type III) and obtain constraints on the matrix elements by experimental data. 
For example, FCNC at the tree level could be compatible with experiments in 
certain regions of the parameter space if we consider the exchange of heavy 
scalar or pseudoscalar Higgs fields |I3 or by cancellation of large contributions 
with opposite sign. 

Finally, a very important test of FC vertices in leptons could be provided by 
the (<? — 2) /2 muon factor (a M ). Recent measurements of a p at BNL [p5| have 
reduced the uncertainty by a factor of two from previous measurements. On the 
other hand, SM estimations on come from QED, electroweak, and hadronic 
contributions pqj; the hadronic ones carry the bulk of the uncertainties. In 



2 



particular, the light by light hadron contributions have been recently corrected 
[p7j. From the experimental and SM values of a M it is possible to estimate 
the room for new physics (Aa M ) available for this observable. Several works 
constraining LFV processes from previous estimations on Aa M have been carried 
out [^8| In order to constrain such kind of processes, we shall use the 

following interval for Aa p at 95% confidence level (C.L.), reported by |3Cj 

9.38 x 1CT 10 < Acifi < 51.28 x 1CT 10 (1) 

In this paper we examine constraints on the Flavor Changing (FC) vertices 
in the lepton Yukawa sector of the 2HDM type III. Such bounds will be gotten 
from Aa^ and the leptonic decays [i — > &y, t — > fij and r — > fififi. In addition, 
we predict upper limits for the decays r — > ej and r — > eee. 



2 LFV Processes 



The Yukawa interactions in the 2HDM type III, relevant for our purposes are 

- £ Y = —^—EMeE (cos aH° - sin ah ) + -^Ef E E (sin aH° + cos ah ) 
2M W v ' V2 K ' 

+ ^ECe^EA +h.c. (2) 

where H°(h°) denote the heaviest (lightest) neutral CP— even scalar, and ^4° is 
a CP— odd scalar. E refers to the three charged leptons E = (e, /i, r) and 
Me, £e are the mass matrix and the matrix of flavor changing vertices respec- 
tively. Finally, a is the mixing angle in the CP— even sector. We are working 
in the simplest parametrization of the 2HDM type III in which only one of the 
doublets acquire a Vacuum Expectation Value (VEV) [jwj . 

Now, we write down the expressions for the processes relevant in our analysis. 
For Aa M we shall use the integral expression for the sake of accuracy. The 
expressions for Aa M , are given by 



Aa M = E^ + E 7 



p 

p 



T r<2 m l f 1 x 2 (l-x±m T /m tl ) 

J S(P) = C S(P)^ / » » , / 2 — -2V-, a,t2 n ~^ dx W 



S{P) 8^ J mix 2 + (m 2 - mj) x + M| (p) (1 - x) 

where Is(p) is an integral involving an Scalar (Pseudoscalar) Higgs boson with 
mass Mg^p), and Cg(p) is the corresponding coefficient in the Yukawa La- 
grangian (||). On the other hand, the expressions for r (t — > /J.J), V (r — * ej) 
and r (fx — ► e7) yield 

r(r->/7j = 4, r — R{m H o,m h o,m A o,a,£ T T) , 



4tt 4 V2 
16tt 



( m 4 m 

r (M->e7) = g r £r r ^rzfa 11 S (m H o,m h o,m A o,a,£, TT ) . (4) 
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where I = e,/i denotes a light charged lepton. In addition, we have defined 

\[2^ TT sin 2 a \ In [m H o / m T ] 



R(m H o,m h o,m A o,a,£ TT ) = 



m-r sin 2a 



m T sin 2a 



'Gf 



1 h° 



\/2£tt cos 2 a \ In [m h o / m T ] 



2 \J~2^ TT In [m^o /7 



'Gf 
i 



'A" 



S (m H o , m h o , m A o , a, £ TT ) 



In [ttiho /m T ] 



sin a * h cos a 



In [m/jo / m T ] 



m 



h° 



In [m^o/m T ] 



A 



(5) 



Finally, the expression for a lepton L going to three leptons of the same flavor 
I is given by 



r (L -> HZ) 



2048tt 3 
6i 



sin 2a 




1 



1 



m, 



■ 2 9 ^ 

sm a cos a 1 



in 



A° 



(6) 



2.1 Obtaining the bounds 

^,From the previous section, the free parameters that we are involved with, 
are: the three neutral Higgs boson masses {m h a , m H o , m A o ) , the mixing angle 
a, and some flavor changing vertices £y. Based on the present bounds from 
LEP2 we shall assume that m^o w 115 GeV. In addition, we shall assume 
that m,A° ^ both assumptions will be held throughout the document. 

Now, since we are going to consider plots in the £y — m A o plane, we should 
manage to use appropiate values of (m H o , a) in order to sweep a wide region of 
parameters. In order to sweep a reasonable set of this couple of parameters we 
utilize for m^o values of the order of 115 GeV (light), 300 GeV (intermediate), 
and very large masses (heavy). As for the angle a, we consider values of a = 
(minimal mixing), a — 7r/4 (intermediate mixing), and a = n/2 (maximal 
mixing). We can check that all possible combinations of (m H o,a) could be 
done by considering five cases 1) when m#o ~ 115 GeV; 2) when m H o ~ 
300 GeV and a = tt/2; 3) when rriffo is very large and a = tt/2; 4) when 
m#° ~ 300 GeV and a — 7r/4; 5) when m^o is very large, and a = n/4. In 
all these cases the mass of the pseudoscalar will be considered in the range 
115 GeV< m A o. 

The first bounds come from Aa M . We use the estimated value of it, given 
by ]|(| at 95% C.L. eq. (|j). Since Aa^ in eq. (|l|) is positive, lower and 
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Bounds on £ 2 r 


Bounds on g T Q r 


Bounds on £ 2 T 


case 1 


7.62 x 1CT 4 < Q T £ 8.31 x 10~ 3 


e eT Q T < 7-33 x 10- iS 


£ 2 T < 4.82 x 10- i& 


case 2 


1.29 x lO" 3 < Q T < 4.42 x 10~ 2 


g T Q T < 2.24 x l(T ib 


( z eT < 2.77 x 10~ i4 


case 3 


1.53 x 10- a < Qr 


g T Q T < 2.24 x 10- ib 


Ct T < 2.76 x 10~ i4 


case 4 


9.57 x 10~ 4 < Clr ^ L4 ° x 


< 2.10 x 10- 1 ' 


& < 8.22 x 10~ ib 


case 5 


1.02 x 10~ 3 < £ 2 T < 1.66 x 10~ 2 


< 2.93 x 10- 1 '' 


l'i r < 9.65 x 10~ ib 



Table 1: Constraints on the mixing parameters £ 2 T , CerCur an(I £er f° r tnc 
five cases mentioned in the text. The two former are generated from Aa M and 
r (/i — * ej) respectively, while the latter comes from the combination of the 
lower limit on £ 2 T and the upper bound on £ 2 T £ 2 T . 

upper bounds for the FC vertex £ MT can be gotten at 95% C.L. The results are 
indicated in table (|l]) column 1. The lower bounds in each case are obtained 
when m^o « 115 GeV, and using the minimum value of Aa M in eq. ([!]), while 
the upper bounds are obtained when ^4° is very heavy and using the upper limit 
for Aa M in eq. ([!]). From these results a quite general and conservative allowed 
interval can be extracted 

7.62 x 10~ 4 < £ 2 r < 4.44 x 10~ 2 . (7) 

Furthermore, upper bounds for the product Q, T ^ T are obtained from the 
expression of the decay width r (/i — > ej) in eq. (0) and from the experimental 
upper limit T (/x — > ej) < 3.6 x 10~ 30 GeV p2]| . The most general upper bounds 
are obtained for A very heavy. The results are shown in table (0) column 2. 
From this table we infer that quite generally, the upper limit is 

e e A% < 2.24 x 10- 16 (8) 

Moreover, combining these upper limits with the lower bounds on £ 2 T given 
in the first column of table (|TJ) , we find upper limits on £ 2 T . The results appear 
on table (Q) third column and the general bound can be written as 

f eT < 2.77 x 10~ 14 . (9) 

Noteworthy, these constraints predict a strong hierarchy between the mixing 
elements £ MT and £ er i.e. |£ e i-| << |£a»t| and they differ by at least five orders of 
magnitude. 

On the other hand, from eq. (|4|) we see that the decay widths r (t — > /17) 
and r (t — * /i/i/i) depend on two mixing vertices £ 2 r , £ rT and £ 2 T , £ MM respec- 
tively. Then, we can get conservative constraints on the diagonal mixing ver- 
tices £ TT , £ MM by using once again the lower bounds on £ 2 T obtained from 
Aa M . Since r (t — > /17), and r (r — * fifi^) are rather complicate functions of 
£ rT , and respectively, we present these constraints in the form of contour- 
plots in the m^o — £ TT plane and the m^o — £ MM plane, figures [l] and S each one 
for the five cases. 
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Figure 1: Contourplots in the £ Tr — m^o plane for each of the five cases cited 
in the text. On left: case 1 corresponds to the long-dashed line, case 4 to the 
short-dashed line, and case 5 to the solid line. On right: case 2 corresponds to 
dashed line, and case 3 is solid line. 




Figure 2: Contourplots in the £ w — m^o plane for each of the five cases cited 
in the text. On left: case 1 corresponds to the long-dashed line, case 4 to the 
short-dashed line, and case 5 to the solid line. On right: case 2 corresponds to 
dashed line, and case 3 is solid line. 



We can check that for each contourplot there is a value of m^o for which 
£ Tr or stays unconstrained, and they are shown in tables |], and ||| respec- 
tively. Additionally, the bounds for £ TT , £ w when m^o is very large and when 
m^o w 115 GeV are also included in tables pi || respectively. We see that the 
general constraints for £ TT read 

- 1.8 x 1(T 2 < £ TT < 2.2 x 1(T 2 , 

-1.0 X 10~ 2 < Crr < 1-0 X 10~ 2 (10) 

for rriA° very large and for rriA° sa 115 GeV respectively. These constraints are 
valid for all cases, except for the third case with to^o very large, since in that 
scenario £ Tr remains unconstrained. Now, for £ w the general bounds read 

IUI £ 0.12, 

IUI £ 0.13 (11) 

for m^o very large and for m^o m 115 GeV respectively. Once again these 
constraints are not valid for the third case when m^o is very large, but are valid 
in all the other cases. 

The presence of diagonal mixing vertices could play a crucial role in looking 
for FCNC. At this respect, it is worthwhile to point out that the relative cou- 
plings (the quotient between the Yukawa couplings of the new physics and the 
SM Yukawa couplings) in model type III 



(fe°7/J y/2M w £f f 

-7 f- = cos a - — 

sm ) gm f 

\ 9 h a lf) , VlMwtff 

-, t- = - sm a + — 



-—sin a + cos a (12) 

sm ) gmf 

are not universal because of the contribution of the factor C/// TO /- As a manner 
of example, the relative couplings H°tt, H /!/!, H°ee are in general different, a 
deviation from the universal behavior of the relative couplings of the type Hff 
could be a clear signature of FCNC at the tree level. In the case of the 2HDM 
type I, all the relative couplings of the form Hff for a certain Higgs are equal. 
In the model type II, relative couplings of this form are equal for all up-type 
fermions, and for all down-type fermions. By contrast, in the model type III all 
relative couplings of the form Hff can be different each other. In addition, in 
the 2HDM with FCNC the pseudoscalar Higgs couples to a pair of fermions of 
the same flavor by means of the matrix element £//, different from the case of 
the 2HDM with no FCNC in which the pseudoscalar couples through the mass 
of the corresponding fermion. 

We can check that our bounds on £ Tr still permits either a huge suppression 
or a huge enhancement on the couplings Htt with H being any neutral Higgs. 
For instance, in the case 5 by using £ rr ~ 0.01 the coupling h°TT vanishes, while 
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Values of m A o for 
£ rT unconstrained 


£ rr intervals for 
m A o very large 


£ TT intervals 

for m A o ~ 115 GeV 


case 1 


m A0 sa 170 GcV 


-0.0072 < Crr < 0.0072 


-0.010 < Crr < 0.010 


case 2 


m A o as 440 GeV 


-0.018 < i TT < 0.018 


-0.0043 < £, TT < 0.0043 


case 3 




unconstrained 


-0.0036 < Crr < 0.0036 


case 4 


m A o « 228 GeV 


-0.0075 < i TT < 0.022 


-0.0094 < i TT < 0.0036 


case 5 


m A o w 250 GcV 


0.00024 < £ TT < 0.021 


-0.0093 < Crr < 0.0026 



Table 2: Bounds extracted from the contourplots shown in fig. |l| The first 
column indicates the values of m A o for which £ TT stays unconstrained for each 
of the five cases. The second and third columns show the allowed intervals on 
£ TT when m A o is very large and when m^o rj 115 GeV respectively. 





Values of m A o for 


^jy, intervals for 


^juju intervals 




£ PM unconstrained 


m^o very large 


for m^o w 115 GeV 


case 1 


m A o as 115 GeV 


IWI < 0.043 


unconstrained 


case 2 


m A o « 300 GeV 




IWI < 0.055 


case 3 




unconstrained 


\^}< 0-043 


case 4 


m A o w 152 GeV 


If^l < 0.058 


<0-13 


case 5 


m^o s» 163 GeV 


1^1 < 0.061 


IWI<o.n 



Table 3: Bounds extracted from the contourplots shown in fig. |[ The first 
column indicates the values of m A o for which £ w stays unconstrained for each 
one of the five cases. The second and third columns show the allowed intervals 
on when m A o is very large and when m A o as 115 GeV respectively. 
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using £ Tr ~ 0.1, we see that the contribution coming from the term proportional 
to £ TT , is about 10 times larger in magnitude than the contribution coming from 
the term proportional to the tau mass. In the case 5 both values £ TT = 0.01 and 
£ rT = 0.1 are allowed at least for m^o around 250 GeV, as it is shown in fig. 
([[]). However, for a very light (m A o w 115 GeV), or a very heavy pseudoscalar 
Higgs boson, constraints on £ rT are considerably stronger in all cases, see table 
||. Similarly, strong or weak couplings of H°tt and/or A°rr are still allowed. 
Of course, the same pattern is accomplished by the bounds for since they 
are weaker than the constraints on £ TT . 

Based on the bounds obtained above, we are able to estimate upper limits for 
some leptonic decays by using the expressions (H), and (||). We shall assume 
that |£ Tr | « |£ ee | < 0.1, from this assumption we can obtain the following upper 
bounds 

r(r^e 7 ) < 1.5 x 10^ 27 , 

r(r^eee) < 5 x 10~ 29 . (13) 

If we compare with the current experimental upper bounds r (r — > ej) < 6.12 x 
10 -18 , r (r — > eee) < 6.57 x 10~ 18 we see that these decays are predicted to 
be very far from the reach of next generation experiments in the context of the 
2HDM type III, unless that £ TT , and/or £ ee acquire unexpectedly large values. 



3 Concluding remarks 

We calculated some constraints on lepton flavor violating vertices in the frame- 
work of the 2HDM type III, by assuming m h o ss H5GeV and m^o > rn h o. Since 
the most recent estimations of Aa^, still provides an important window for new 
Physics, we obtain from it lower and upper bounds on the mixing vertex ^ r at 
95% C.L. Specifically, an allowed interval of 7.62 x 10~ 4 < £ 2 T < 4.44 x 10~ 2 was 
found in a quite wide region of parameters. Of course, we should realize that 
both SM test and experimental measurements of are still being scrutinized 
and current results are not definitive at all. However, if not severe changes 
occur in forthcoming experiments and SM estimations, these constraints could 
continue being valid at least at a lower confidence level. Future improvements 
on both estimations should elucidate this point. 

Based on these constraints, and on the leptonic decays \i — > &y, t — > hj, and 
r — > fififi we got the following conservative and quite general bounds on LFV 
vertices 

i\ T < 2.77 x 10~ 14 . (14) 

For £ TT we find 

-1.8xl0~ 2 < £tt < 2.2 x 10~ 2 , 

|£r T | < 1.0 X 10" 2 (15) 
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for niAo very large and for m^o w 115 GeV respectively. For we find 

IUI £ °- 15 ( 16 ) 

for m^o very large and for m^o w 115 GeV respectively. Notwithstanding, those 
bounds on £ Tr and £ w are considerably relaxed for certain specific values of 
m^o. In that case, the room available for them is such that they permit either 
a strong enhancement or a strong suppression of the couplings Htt and/or 
HfiJI for any neutral Higgs boson. Additionally, according to these constraints 
we find that |£^ T | is at least five orders of magnitude larger than |£ eT | revealing 
a strong hierarchy between them. 

Furthermore, we estimate upper limits on the decay widths T(t — ► ey), 
and r(r — ► eee), finding that they are basically hopeless to look for LFV pro- 
cesses in near future experiments, at least in the framework of the 2HDM type 
III with heavy Higgs bosons. 
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